Algorytmy danych geoprzestrzennych

Obliczenia geodezyjne

Krzysztof Dyba

Jaki kształt ma Ziemia?

Blue Marble 2012, Suomi NPP

Geoida

International Centre for Global Earth Models

https://essd.copernicus.org/articles/11/647/2019/

Geoida

Geoida to model kształtu Ziemi, który przedstawia złożoną, nieregularną powierzchnię, na którą siła grawitacji działa wszędzie prostopadle do powierzchni.

Geoida nie ma matematycznie zdefiniowanego kształtu. Jej kształt wynika z różnic w ziemskim polu grawitacyjnym spowodowanymi zmianami w gęstości i rozkładzie masy Ziemi (np. występowanie gór i dolin).

Źródłem danych zazwyczaj są satelitarne pomiary grawitacyjne (np. misja GRACE).

Matematyczne reprezentacje

Sfera (góra)

Elipsoida obrotowa (lewo dół)

Elipsoida (prawo dół)

Wikipedia

Elipsoida

Elipsoida to obiekt geometryczny, w którym wszystkie płaskie przekroje są elipsami lub okręgami. Wszystkie trzy osie mogą mieć różne długości.

\[ \frac{x^2}{a^2} + \frac{y^2}{a^2} + \frac{z^2}{b^2} = 1 \]

Elipsoida obrotowa (sferoida)

Elipsoida obrotowa to obiekt geometryczny utworzony przez obrót elipsy wokół jednej z jej osi. Jeśli elipsa zostanie obrócona wokół mniejszej osi, to w rezultacie powstanie elipsoida spłaszczona, która przypomina kształem Ziemię.

Jest to szczególny przypadek elipsoidy, w której dwie z trzech osi są równe.

\[ \frac{x^2 + y^2}{a^2} + \frac{z^2}{c^2} = 1 \]

Jeśli:

\(c < a\) – elipsoida spłaszczona

\(c > a\) – elipsoida wydłużona

Sfera

Sfera to idealnie okrągły obiekt geometryczny, w którym każdy punkt na powierzchni jest w równej odległości od środka. Jest to podstawowe i ogólne przybliżenie kształtu Ziemi, ponieważ nie uwzględnia spłaszczenia na biegunach.

Jest to szczególny przypadek elipsoidy, w której wszystkie trzy osie mają równą długość.

\[ x^2 + y^2 + z^2 = r^2 \]

Obrazowanie Ziemi na mapie

Podstawowy problem

Przedstawienie Ziemi na mapie obejmuje rzutowanie trójwymiarowego kształtu na dwuwymiarową powierzchnię. Ponieważ powierzchnia Ziemi jest zaokrąglona, proces ten wprowadza zniekształcenia związane ze kształtem, powierzchnią, odległością czy kierunkiem.



Jakie jest rozwiązanie?

Odwzorowanie kartograficzne

Odwzorowanie kartograficzne to matematyczna transformacja szerokości i długości geograficznej z powierzchni sfery lub elipsoidy na płaszczyznę.

Podstawowe typy odwzorowań:

  • Walcowe – Ziemia rzutowana jest na cylinder, który następnie jest rozwijany tworząc płaską mapę.
  • Stożkowe – Ziemia rzutowana jest na stożek, który jest rozwijany.
  • Azymutalne (płaszczyznowe) – Ziemia rzutowana jest na płaską płaszczyznę, wyśrodkowaną w określonym punkcie.

Przykład

Przykład

Elipsy Tissota

Wnioski

  • Przedstawianie Ziemi na mapie to skomplikowany proces wymagający kompromisów.
  • Każde odwzorowanie ma swój własny zestaw zniekształceń.
  • Różne cele wymagają odmiennych podejść.
  • Odwzorowanie jest zawsze przybliżeniem.

Obliczenia

Odległość

Odległość euklidesowa: 8116 km
(Odwzorowanie azymutalne równoodległościowe)

Odległość sferyczna: 6855 km

Odległość euklidesowa

Odległość euklidesowa to miara odległości w linii prostej między dwoma punktami obliczana z następującego równania:

\[ d = \sqrt{(\Delta x)^2+(\Delta y)^2} = \sqrt{(x_2 - x_1)^2 + (y_2 - y_1)^2} \]

Wynika ona bezpośrednio z twierdzenia Pitagorasa o trójkątach prostokątnych:

\[ a^2 + b^2 = c^2 \]

gdzie \(a\) i \(b\) to długość przyprostokątnych, a \(c\) to długość przeciwprostokątnej.

Odległość sferyczna

Odległość sferyczna (ortodroma) to najkrótsza odległość między dwoma punktami na powierzchni sfery.

Wzory:

  1. Równanie haversine
  2. Z długości cięciwy

Odległość elipsoidalna

Odległość powierzchniowa (z NMT)

Odległość powierzchniowa (z NMT)

Powierzchnia

Rzeczywista wielkość komórek

Precyzja współrzędnych

https://xkcd.com/2170